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摘 要 : 为 实现 定量 化 分 析 温 、 光 和 水 资源 对 果树 生长 的 影响 ， 本 研究 以 成 龄 骏 束 树 为 研究 对 象 ， 提 出 了 


基于 校正 WOFOST 模 型 的 束 树 生长 和 水 分 运 移 模拟 方法 。 利 用 2016 和 2017 年 的 田间 试验 观测 数据 ， 重 点 校 
正 WOFOST 模 型 的 物候 学 发 育 、 初 始 化 、 绿 叶 、C0, 同 化 、 干 物质 分 配 、 呼 吸 作用 和 水 分 利用 参数 。 在 田间 


尺度 ， 完 成 总 地 上 生物 量 (TAGP)、 叶 面积 指数 (LAI) 和 土壤 水 分 含量 的 动态 模拟 和 精度 验证 ; 在 县 域 尺 
E, 使 用 55 个 果园 的 最 大 LAI1、 单 产 、 实 际 燕 散 量 (ET,) 和 水 分 利用 效率 (WUE) 数据 评价 模型 区 域 尺 度 
的 模拟 性 能 。 结 果 表 明 ， 在 田间 尺度 ,校正 模型 模拟 不 同 灌溉 梯度 TAGP 的 决定 系数 尼 范 围 为 0.92 一 0.98 ， 


响 的 定量 化 分 析 提 供 新 思路 。 


归 一 化 均 方 根 误差 NRMSE H 8.7% ~20.5% ; PEF LATAS] R^ 35 EIH 0.79~0.97, NRMSE J 8.3%~21.1%; 模拟 
土壤 水 分 含量 的 决定 系数 尼 范 围 为 0.29~-0.75，NRMSE 为 4.1%~-6.1%。 在 县 域 尺度 ， 两 年 模拟 最 大 LAI 与 实 
测 LAI 的 尼 分 别 为 0.64 和 0.78，NRMSE 分 别 为 13.3% 和 10.7%; 模拟 单产 的 尼 分 别 为 0.48 和 0.60，NRMSE 分 
HE 12.1% 和 11.9%; 模拟 ET, 均 方 根 误差 分 别 为 36.1 mm (7.9%) 和 30.8 mm (7.4%); 模型 也 表现 了 较 高 
的 WUE 模 拟 精度 (10%<NRMSE<20%)， 均 方 根 误差 RMSE 值 分 别 为 0.23 和 0.28 kg/m*。WOFOST 模 型 在 田间 
和 县 域 尺度 都 取得 了 较 高 的 束 树 生长 和 水 分 运 移 模 拟 精度 ， 可 为 土壤 、 气 象 、 灌 溉 管理 和 讽 树 生长 耦合 影 
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1 引言 


ARAR EAE, BER. ERIT VICE 
营养 成 分 ， 具 有 重要 的 食用 和 医学 价值 。 中 国 红 
束 种 植 区 域 主要 有 新 疆 、 山 西 、 陕 西 、 河 北 、 河 
南 和 山东 等 地 区 。 新 疆 南 部 地 区 光 热 资源 丰富 、 
日 照 时 间 长 、 昼 夜 温 差 大 ， 是 中 国 最 重要 的 优质 


收 稿 日 期 2021-03-23 ”修订 上 日期: 2021-04-20 


红 束 生产 基地 ， 年 产量 372.8 万 吨 ， 占 全 国 总 产 
量 的 50% ,已 成 为 南 疆 农民 增收 和 人 口 肾 集 的 
重要 支柱 产业 。 然 而 ， 新 疆 地 区 夏季 温度 高 、 落 
发 蒸腾 量 大 、 王 旱 少雨 ， 水 资源 紧缺 ， 如 何 定量 
评价 蒸 散 和 水 分 利用 效率 是 干旱 区 农业 节 水 的 关 
键 技术 问题 。 现 有 关于 束 树 生长 水 效应 和 灌溉 管 
理 研 究 主要 集中 在 灌水 处 理 对 产量 和 品质 的 影 
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响 中 、 管 理 制度 对 土壤 水 分 的 影响 ”、 灌 溉 模式 
优化 水 肥 管理 ”、 灌 溉 和 施肥 对 产量 的 耦合 影 
响 “等 。 然 而， 这 些 研究 大 多 基于 田间 试验 数 
据 ， 灌 溉 方案 通常 只 适用 于 特定 土壤 和 气候 环 
境 ， 无 法 实现 定量 化 分 析 灌 溉 方案 对 红 束 产量 和 
品质 的 影响 。 

作物 生长 模型 是 基于 生理 和 生态 学 原理 ， 并 
使 用 数学 模型 量化 受气 候 、 土 壤 和 管理 条 件 影响 
的 作物 生长 发 育 过 程 。 在 过 去 的 儿 十 年 中 , 已 经 
开发 的 几 个 比较 成 熟 的 作物 模型 包括 WOFOST 
(WOrld FOod Studies) 、DSSAT (Decision Sup- 


port System for Agrotechnology Transfer) #1] EPIC 
(Environmental Policy Integrated Climate) F *, 
这 些 模型 各 具 特 点 和 优势 ， 但 通常 都 包含 了 物候 
学 、CO, 同 化 (光合 生产 )、 呼 吸 作 用 、 同 化 物 
分 配 和 土壤 水 循环 等 模块 ， 被 广泛 应 用 于 农业 生 
产 管理 以 及 与 社会 相关 的 应 用 中 ， 如 气候 变化 响 
应 ”， 作 物产 量 预 测 ”以 及 环境 、 水 分 和 作物 
生长 的 耦合 作用 等 。 已 有 的 作物 生长 模型 主 
要 聚焦 一 年 生 的 大 田 作 物 ， 关 于 果树 作物 模型 研 
究 方面 ，Fishman 和 Geénard "基于 水 和 干 物质 
运输 机 制 开发 了 桃 树 果 实生 长 模型 计算 干 旱 和 
新 增 质量 变化 所 表示 的 季节 性 和 昼夜 水 果 生 长 情 
况 ， 以 分 析 作物 负荷 下 的 水 分 胁迫 条 件 。Lopez 
4B 8) 基于 形态 学 开发 了 桃 树 生长 的 计算 机 模拟 
模型 ， 实 现 桃 树 生 长 的 三 维 描述 ， 可 用 于 树木 修 
剪 管理 。Lescourret 等 ”用 全 局 模型 解释 了 上 肉 性 
和 雄性 蕨 的 开花 、 授 粉 和 果实 生长 之 间 的 关系 ， 
并 集成 了 生产 管理 输入 ,仿真 了 果园 管理 对 猕猴 
桃 果 实数 量 和 尺寸 的 影响 ， 该 模型 对 气候 和 技术 
操作 的 变化 都 很 敏感 。Green 等 "开发 了 果树 冠 
层 蒸 腾 和 光 能 截获 模型 ， 用 于 证 明 蔡 腾 作 用 如 何 
响应 气象 环境 的 变化 ， 敏 感性 分 析 表 明 ， 光 拦截 
受 叶 面积 和 叶 光 学 特性 变化 的 影响 最 大 ， 而 蒸腾 
作用 受 叶 面积 和 叶 电 导 变 化 的 影响 最 大 。Costes 
等 "使 用 混合 随机 和 生物 力学 模型 模拟 了 果树 
发 育 。 

然而 ， 已 报道 的 果树 生长 模型 缺少 生理 生化 


过 程 的 机 理性 描述 。WOFOST 是 一 种 机 理 模型 ， 
可 以 根据 物候 学 、CO, 同 化 、 蒸 腾 作 用 、 呼 吸 作 
用 以 及 这 些 过 程 如 何 受 环境 条 件 影 响 来 解释 作物 
生长 “"。 该 模型 的 过 程 描述 是 通用 的 ， 通 过 参数 
校正 可 用 于 不 同 作 物 的 生长 模拟 "7" 和 水 分 限制 
产量 评估 ,不 仅 被 广泛 用 于 定量 分 析 一 年 生 
大 田 作物 的 生长 过 程 ,而 且 具 有 束 树 生长 模拟 
和 产量 评估 的 潜力 。 基 于 WOFOST 实 现 多 年 
生 果 树 作物 的 水 分 限制 生长 模拟 ， 进 而 实现 土 
培 、 气 象 、 灌 溉 管理 和 果树 作物 生长 耦合 影响 的 
定量 化 分 析 ， 有 望 使 果树 灌溉 管理 从 试验 统计 分 
析 匹 入 定量 、 机 理性 描述 阶段 。 因 此 ， 本 研究 以 
矮 化 密植 骏 束 为 研究 对 象 ， 在 WOFOST 模 型 作 
物 参数 校正 基础 上 ， 实 现 吏 树 生 长 和 水 分 运 移 模 
拟 过 程 。 并 使 用 大 田 试 验 测量 的 总 地 上 生物 量 
(Total Above-Ground Production，TAGP) 、 叶 面 
积 指 数 (Leaf Area Index, LAI) 和 土壤 水 分 含 
量 验证 田间 尺度 的 模拟 精度 ; 使 用 县 域 尺 度 测 量 
的 55 个 果园 的 最 大 叶 面 积 指数 、 单 产 、 实 际 薰 
散 量 和 水 分 利用 效率 评价 校正 WOFOST 模型 在 
县 域 尺 度 不 同 束 园 的 模拟 性 能 。 


2 材料 与 方法 


2.1 研究 区 概况 


研究 区 位 于 新 疆 维吾尔 自治 区 阿拉 尔 市 
(80930'E—81958'E, 40922'N—40957'N), ZL fi 
植 面 积 约 4.5 万 公 质 ， 约 占 全 国 红枣 产 量 的 18. 
属于 暖 温 带 极端 大 陆 性 干旱 荒漠 气候 ，2016 年 和 
2017 年 平均 降雨 量 分 别 为 106 和 98 mm， 主 要 分 
布 在 夏季 ， 太 阳 辐 射 年 均 133.7—146.3 kcal/em’, 
4E. H HE 2556.3—2991.8 Ph， 年 平均 气温 在 10.8~ 
12.5 "C 之 间 ， 最 大 日 温差 为 20 *C。 丰 富 的 光 热 
资源 为 优质 红枣 种 植 提供 了 良好 的 自然 条 件 。 
2.2 试验 设计 

在 田间 试验 区 域 (图 1 蓝 色 方块 定位 区 域 ) 


选择 成 龄 束 园 在 2016 和 2017 年 开展 灌溉 试验 ， 
采用 滴灌 方式 在 生育 期 设计 $ 个 灌溉 梯度 (不 包 
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括 冬 灌 和 春灌 ): 灌溉 定额 CK (375mm, D1), 
90% CK (338 mm, D2), 80% CK (300 mm, 
D3), 70% CK (263 mm, D4) 和 60% CK 
(225 mm, D5). iE WEE 812 5 T 1H Hh ee a t BR 
量 ， 生 长 季 灌 溉 10 次 ， 各 小 区 施肥 量 和 其 他 农 
艺 措施 相同 。 
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图 1 试验 和 采样 表 园 分 布 


Fig. 1 Distribution of experimental and sampled jujube gardens 


2.3 测试 项 目 及 方法 


在 田间 试验 区 域 测量 的 主要 数据 包括 以 下 
内 容 。 

D 物候 学 发 育 时 间 : 记录 萌芽 ， FE (e 
果 ) 和 成 熟 日 期 。 

@) 干 重 测量 : 每 10 d 取 样 一 次 ， 生 长 季 取 
样 10 次 ， 每 次 选 均 匀 的 束 树 3 棵 ， 取 全 树叶 片 、 
新 枝 和 果实 带 回 实验 室 在 80 °C 下 烘 干 至 恒 重 后 
测量 各 器 官 干 重 。 

(3) 光合 有 效 辐射 和 叶 面 积 指数 : 每 10 d 使 
用 SunScan 冠 层 分 析 仪 (Delta 公司， 英国 ) 测定 
不 同 冠 层 空间 的 有 效 叶 面 积 指数 。 

@ 光合 作用 : 每 10 4 测量 一 次 ,使 用 
LI6400XT 便携 式 光 合作 用 测试 仪 分 层 测试 净 光 
合 速率 、 气 孔 导 度 、 胞 间 CO, YR EE AZ ES HER 

O 土壤 水 分 : 灌溉 前 后 在 0~100 cm 土壤 深 
度 每 20 cm 用 根 钻 取 土 样 ， 测 定 土壤 容积 含水 
量 ， 每 年 测量 14 次 。 

© 其 他 土壤 参数 : 土壤 田 持 、 容 积 密度 、 
饱和 土壤 含水 率 、 土 壤 水 响应 曲线 和 渗透 系数 等 
直接 取样 带 回 实验 室 进行 测量 ,同步 记录 灌溉 日 
期 和 灌溉 量 。 


QD 气象 数据 : 由 塔里木 大 学 红枣 研究 基地 
架设 的 小 型 气象 站 提供 ,包括 15 min 间 隔 的 温 
度 、 湿 度 、 辐 射 、 风 速 、 降 雨 和 气压 数据 等 。 

在 县 域 尺度 观测 55 个 束 园 (图 1 红 点 坐标 ) 
的 初始 束 树 干 重 、 最 大 发 育 阶段 的 叶 面 积 指数 、 
最 终 产 量 和 实际 蒸 散 。 叶 面积 指数 观测 时 间 分 别 
为 2016 年 7 月 24 日 和 2017 年 7 月 27 日 ， 每 年 11 
月 份 观测 束 园 产量 。 实 际 蒸 散 量 ET, 计算 如 公 
式 Cl) Bs", 

ET,=P,+I- D -ASM (1) 

HE, POL D 和 ASM 分 别 为 有 效 降 雨 
量 、 灌 溉 量 、 深 层 土壤 渗透 、 生 长 季 开 始 和 结束 
时 土壤 平均 水 分 含量 差 值 。 新 疆 地 区 骏 束 树 根 系 
存在 于 0~100 cm 7". (主要 分 布 在 垂直 60 cm 内 
XE), ， 研 究 中 测量 100 cm 土壤 层 渗透 量 作为 深 
层 土壤 渗透 量 。 研 究 区 地 下 水 较 深 ， 地 表 径 流量 
极 少 ， 所 以 研究 中 未 考虑 地 下 水 和 地 表 径 流 对 
ET7, 的 影响 。 

水 分 利用 效率 (Water Use Efficiency, 
WUE) 是 节 水 农业 的 重要 指标 ， 尤 其 在 干旱 和 
半 干 旱地 区 ， 本 研究 中 也 将 使 用 WUE 进一步 验 
证 校正 模型 的 区 域 尺 度 模拟 精度 。WUE 等 于 实 
测 产 量 与 57, 的 比值 。 


2.4 WOFOST 模 型 参数 校正 


WOFOST 模 型 参数 包括 气象 、 土 壤 和 作物 
参数 ， 由 于 研究 区 域 为 县 域 尺度 ,气象 数据 空间 
变异 性 较 小 ， 因 此 直接 使 用 气象 站 观测 数据 ， 土 
壤 基 本 属性 参数 直接 使 用 测量 值 。 本 研究 主要 校 
正 作物 参数 ， 对 于 敏感 性 较 低 的 参数 ， 如 发 芽 的 
最 低温 度 TBASEM 和 叶片 老化 的 最 低温 度 
TBASE 通过 查阅 文献 “确定 ， 敏 感性 较 高 的 参 
数 ， 如 初始 生物 量 TDWI 和 叶片 死亡 率 SPAN 通 
过 2016 年 D1 灌溉 处 理 的 田间 观测 数据 进行 校 
正 。 参 考 Wit 等 2 的 关于 WOFOST 模 型 的 校准 
方法 ,结合 束 树 生理 特点 ，WOFOST 模 型 主要 
参数 的 校正 过 程 如 下 。 

(1) 物候 学 参数 校正 。 模 型 校准 从 物候 学 参 
数 开始 ， 通 过 有 效 积温 法 确定 所 。 吏 树 开 始 发 育 
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的 最 低 阅 值 温度 (Lower Threshold Temperature 
for Emergence, TBASEM) 被 设置 为 10 °C 7", 
WOFOST 给 出 的 萌芽 时 最 高 有 效 温 度 (Max Ef- 
fective Temperature for Emergence, TEFFMX) 范 
EIX 18~32°C. HH 2016 Fil 2017 FERH ZH] 
最 高 温度 低 于 30 "C， 萌 芽 期 并 无 高 温 胁迫 ， 
此 ， 将 最 高 有 效 温 度 设置 为 30"C。 从 开始 发 育 
至 萌芽 的 有 效 积 温和 (Temperature Sum from 
Sowing to Emergence, TSUMEM)., WF EFH 
的 有 效 积 温和 (Temperature Sum from Emer- 
gence to Anthesis, TSUMI) 以 及 开花 至 成 熟 的 
有 效 积 温和 (Temperature Sum from Anthesis to 
Maturity，TSUM2) 3 个 指标 通过 观测 每 日 平均 
温度 和 萌芽 、 开 花 以 及 成 熟 日 期 进行 校准 。 有 效 
积温 的 每 日 增加 量 (Daily Increase in Tempera- 
ture Sum, DTSMTB) 根据 最 高 适宜 温度 
(35.5 °C) ”和 最 低 发 育 温度 (10 °C) 7 计算 。 

(2) 初始 干 重 校正 。 现 树 与 一 年 生 作 物 不 
同 ， 如 果 考 虑 上 一 年 的 荃 ， 可 能 会 导致 初始 作物 
干 重 值 过 高 。 本 研究 中 初始 束 树 干 重 被 定义 为 初 
始 新 器 官 (初始 芽 和 根 ) 的 重量 ， 通 过 模拟 和 实 
测 的 总 地 上 生物 量 和 叶 面 积 指数 值 进行 初始 干 重 
TDWI 的 校正 。 县 域 尺 度 的 55 个 果园 的 初始 束 树 
干 重 值 根据 实测 的 芽 重 和 芽 的 分 配 系 数 计算 获 
得 。 其 中 芽 的 分 配 系 数 来 自 于 田间 试验 的 测量 
数据 。 

(3) 绿叶 参数 校正 。 在 所 有 试验 小 区 中 ， 萌 
芽 时 〈 定 义 为 芽 上 第 五 片 叶子 展开 时 ) 的 叶 面 积 
指数 (Leaf Area Index at Emergence, LAIEM) 
的 测量 值 均 小 于 0.004。 另 外 ， 仿 真 结 果 表 明 ， 
当 LAIEM 从 0.0007 (WOFOST 给 出 的 最 小 值 ) 
增加 到 0.01 时 ， 总 地 上 生物 量 、 活 时 干 重 、 活 茎 
干 重 、 果 实干 重 和 LAI 的 模拟 结果 差异 极 小 ， 


乎 不 变 ， 因 此 RGRLAI 被 设置 为 最 小 值 0.05。 叶 
面积 与 干 重 比值 (Specific Leaf Area, SLATB) 
根据 测量 的 总 地 上 生物 量 和 叶 面 积 指数 进行 校 
准 。 叶 片 老化 速率 (Life Span of Leaves Growing 
at 35 Celsius，SPAN) 根据 生长 季 后 期 田间 测量 
的 LAI 进 行 了 校正 。 

(4) CO, 同 化 参数 校准 。 分 别 在 束 树 最 佳 发 
育 下 限 温 度 19.5°*C 和 上 限 温度 35.5°C 拟 合 光 响 
应 曲线 ， 拟 合 获 得 最 大 的 CO, 同 化 速率 (Maxi- 
mum Leaf CO, Assimilation Rate, AMAXTB) 和 
光 能 利用 效率 (Light-use Efficiency for Single 
Leaf, EFFTB) 初始 值 ， 再 通过 实测 总 地 上 生物 
量 和 叶 面 积 指 数 ， 对 消光 系数 (Extinction Coef- 
ficient for Diffuse Visible Light， KDIFTB) 、 
AMAXTB 和 EFFTB 进行 了 校准 。 

(5) 干 物质 分 配 参 数 。FSTB (分 配 到 茎 的 
系数 )、FOTB (分 配 到 有 果实 的 系数 ) 和 FLTB 
(分 配 到 叶 的 系数 ) 为 发 育 阶段 的 函数 ， 使 用 测 
量 和 模拟 的 地 上 各 器 官 干 物质 重量 进行 校准 。 

(6) 水 分 利用 参数 校正 。 使 用 测量 的 土壤 水 
分 含量 校正 蒸腾 速率 修正 系数 (Correction Fac- 
tor Transpiration Rate，CFET) 、 水 分 胁迫 敏感 性 
校正 系数 (Correction Factor for Crop Water 
Stress Sensitivity, DEPNR) 、 初 始 根 深 (Initial 
Rooting Depth，RDI) 、 每 日 根深 最 大 增 量 
(Maximum Daily Increase in Rooting Depth, RRI) 
和 最 大 根深 (Maximum Rooting Depth, RD- 
MCR). 


2.5 模型 验证 


模型 验证 包括 两 部 分 。 在 田间 试验 区 域 ， 使 
用 2017 年 D1~D5 处 理 测 量 的 时 间 序 列 的 总 地 上 
生物 量 、LAI 和 土壤 水 分 含量 验证 束 树 生长 和 土 


此 ，LAIEM 参 数 可 以 设置 为 0.0007 至 0.01 之 间 
的 任何 值 。 测 量 的 LAI 最 大 相对 增加 量 (Maxi- 
mum Relative Increase in LAI, RGRLAI) 均 小 于 
0.05， 而 且 当 RGRLAI 从 0.05 增 加 到 默认 最 大 值 
0.5 时 ， 总 地 上 生物 量 和 叶 面 积 指数 模拟 结果 几 


壤 水 分 运 移 的 模拟 精度 。 县 域 测量 的 55 个 果园 
最 大 发 育 阶 段 的 LAI、 单 产 、 实 际 蒸 散 量 和 水 分 
利用 效率 用 于 验证 模型 区 域 模拟 性 能 。 使 用 决定 
系数 (Coefficient of Determination, R°), HIJAR 
误差 (Root Mean Square Error, RMSE) 和 归 一 
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化 均 方 根 误差 (Normalized Root Mean Square Er- 
ror, NRMSE) 评估 田间 试验 处 理 的 模拟 性 能 ， 
{i FAR’, RMSE, NRMSE, #1 Xt fW (Relative NRMSE = (4) 
Bias Error, RBE) 频率 分 布 和 变异 系数 (Coeffi- _ 
cient of Variation, CV) 评价 县 域 尺度 束 园 的 最 RBE(%) = x Yi x 100% (5) 
大 LAI 和 单产 模拟 性 能 。 使 用 RBE 频率 分 布 、 AM 
RMSE Fil NRMSE VE 4fr SBR 2 AM AIK A} Fl FH ARE nee (6) 


的 一 致 性 ，RAMSE 和 NRMSE 表达 测量 值 和 模拟 个 样本 测量 值 ; 7 表示 测量 样本 平均 值 ; 了 表示 
值 相对 误差 和 绝对 误差 ，NRMSE < 10% 表示 极 
高 精度 ，10% < NRMSE < 20% 表示 高 精度 ， 
20% < NRMSE < 30% 表示 中 等 精度 ，NRMSE > 
30% 表 示 低 精度 5。CF 反 映 了 模拟 结果 空间 变 N ean 
异性 ,Cy <4 RRBESH, mecre Te keba BORITEDPURRE, mimi Sees 
100% 表示 中 等 变异 性 ，CV > 100% 表 示 强 变异 es 

性 "中 。RBE 频 率 分 布 表示 样本 被 高 估 和 低估 的 ” 3 结果 与 分 析 

比例 。R*、RMSE、NRMSE、RBE 和 CV 计算 见 公 

st (2) ~ (6). 3.1 模型 参数 校正 


根据 上 述 模 型 参数 校正 方法 ， 校 正 的 束 树 作 


模拟 样本 平均 值 ，SD 表示 模拟 结果 标准 差 ; n 为 
样本 总 数 ， Es 指标 : 总 地 上 生物 量 ， kg/ha; 
叶 面 积 指数 ，mzmm; 土壤 体积 含水 量 ，cmzycnnm ; 


(2) 。 物 参 数 结果 如 表 1 所 示 。 参 数 分 别 来 自 测 量 
fi (m)、 测 量 值 基础 上 的 校正 值 (m-c)、 佑 测 
f (e). BEE (c) 和 文献 参考 值 。 

(3) 


表 1 模型 主要 作物 参数 校正 值 


Table 1 Calibrated values of main crop parameters of the model 


类 别 参数 名 称 描述 值 单位 来 源 
TBASEM 发 芽 的 最 低温 度 10 °C [22] 
萌芽 参数 TEFFMX 发 芽 的 最 高 温度 30 °C e 
TSUMEM 从 开始 到 发 芽 的 有 效 积温 230 °C m-c 
TSUMI 发 芽 到 开花 的 有 效 积温 967 *C/d m-c 
物候 学 参数 TSUM2 d dn uu 960 *C/d m-c 
DTSMTB100 平均 温度 10 °C 时 每 日 有 效 积温 0 *C/d [22] 
DTSMTB355 平均 温度 35.5 °C 时 每 日 有 效 积温 25.5 °C/d [24] 
TDWI 萌芽 时 初始 作物 总 干 重 5.69—25.9 kg/ha m 
初始 参数 LAIEM 萌芽 时 叶 面 积 指数 LAI 0.004 一 一 m 
RGRLAI LAI 每 日 最 大 增加 量 0.05 a! m 
SLATBO000 DVS = 0 时 叶 面 积 与 干 重 比值 0.00165 ha/kg m-c 
绿叶 参数 SLATB55 DVS = 0.55 时 叶 面 积 与 干 重 比值 0.0013 ha/kg m-c 
SLATB100 DVS = 1 时 叶 面 积 与 干 重 比值 0.0013 ha/kg m-c 
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续 表 1 
类 别 参数 名 称 描述 值 单位 来 源 
SLATB200 DVS =2 时 叶 面 积 与 干 重 比值 0.0014 ha/kg m-c 
SPAN 35.5°C 时 叶片 存活 天 数 40~60 d c 
TBASE 叶片 老化 的 最 低温 度 10 $c [24] 
KDIFTBOO DVS =0 时 可 见 光 漫 反射 的 消光 系数 0.8 = m-c 
KDIFTB200 DVS =2 时 可 见 光 漫 反射 的 消光 系数 0.8 = m-c 
EFFTB19.5 19.5°C 时 叶片 光合 作用 效率 0.495 kg-m?-s/(ha-hr-J) m-c 
CO, 同 化 参数 EFFTB355 35.5°C 时 叶片 光合 作用 效率 0.495 kg:m^'s/Cha:hr-J) m-c 
AMAXTB00 DVS =0 时 最 大 CO, 同 化 速率 39.0 kg/(ha-hr) m-c 
AMAXTBI70 | DVS =1.7 时 最 大 CO, 同 化 速率 39.0 kg/(ha-hr) m-c 
AMAXTB200 DVS =2 时 最 大 CO, 同 化 速率 20.0 kg/(ha:hr) m-c 
CVL 同化 物 转化 到 叶 的 效率 0.732 一 一 c 
T 同化 物 转化 为 生物 量 的 CVO 同化 物 转化 到 果实 的 效率 0.780 —— c 
> BOR CVR 同化 物 转化 到 根 的 效率 0.690 — c 
= CVS 同化 物 转化 到 茎 的 效率 0.751 一 一 c 
N Q10 温度 每 升 高 10 °C 呼吸 速率 相对 增长 率 2 kg: CH,O/(kg:d) m 
e RML 十 的 相对 呼吸 速率 0.03 kg-CH,O/(kg-d) m 
© 维持 呼吸 RMO 果实 的 相对 呼吸 速率 0.01 kg:CH,O/(kg:d) m 
CN RMR 根 的 相对 呼吸 速率 0.01 kg-CH,O/(kg:d) a 
e RMS 茎 的 相对 呼吸 速率 0.015 kg-CH,O/(kg-d) m 
ES FLTB00 DVS = 0 叶 分 配 系数 0.67 一 一 m-c 
© FLTBO012 DVS = 0.12 叶 分 配 系数 0.31 — m-c 
CN FLTB022 DVS = 0.22 叶 分 配 系数 0.41 一 一 m-c 
> FLTB032 DVS = 0.32 叶 分 配 系 数 0.55 一 一 m-c 
x FLTB051 DVS = 0.51 叶 分 配 系数 0.40 一 一 m-c 
人 FLTB097 DVS = 0.97 叶 分 配 系 数 0.15 一 一 m-c 
c FLTB100 DVS = 1.00 叶 分 配 系数 0.1 == m-c 
T FLTB145 DVS = 1.45 叶 分 配 系数 0 —— m-c 
(S) FLTB200 DVS = 2.00 叶 分 配 系 数 0 = m-c 
干 物质 分 配 系数 FSTB00 DVS = 0 茎 分 配 系数 0.33 — m-c 
FSTB012 DVS = 0.12 254) fi R% 0.69 一 一 m-c 
FSTB022 DVS = 0.22 茎 分 配 系数 0.59 一 一 m-c 
FSTB032 DVS = 0.32 24} HUA Re 0.45 — m-c 
FSTB051 DVS = 0.51 茎 分 配 系数 0.6 —— m-c 
FSTB097 DVS = 0.97 茎 分 配 系 数 0.85 一 一 m-c 
FSTB100 DVS = 1.00 茎 分 配 系数 0.43 MES m-c 
FSTB145 DVS = 1.45 茎 分 配 系 数 0.2 —— m-c 
FSTB200 DVS = 2.00 葵 分 配 系 数 0 一 一 m-c 
FOTB00 DVS = 0 果实 分 配 系数 0 一 一 m-c 
FOTB012 DVS = 0.12 果 实 分 配 系数 0 —— m-c 
干 物质 分 配 系数 FOTB022 DVS = 0.22 果 实 分 配 系 数 0 一 一 m-c 
FOTB032 DVS = 0.32 果 实 分 配 系 数 0 一 一 m-c 
FOTBO051 DVS = 0.51 果 实 分 配 系数 0 一 一 m-c 
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ES 参数 名 称 描述 值 单位 来 源 
FOTB097 DVS = 0.97 果 实 分 配 系数 0 一 一 m-c 
FOTBI100 DVS = 1.00 果 实 分 配 系 数 0.47 一 一 m-c 
FOTB145 DVS = 1.45 果 实 分 配 系 数 0.8 一 一 m-c 
FOTB200 DVS = 2.00 果 实 分 配 系数 1 一 一 m-c 
RDRSTBO00 DVS = 0 茎 相对 死亡 率 0 一 一 e 
死亡 率 
RDRSTB200 DVS =2.0 葵 相对 死亡 率 0 一 一 e 
CFET 蒸腾 速率 校正 系数 1.02 一 一 c 
DEPNR 作物 水 分 胁迫 敏感 性 校正 系数 1.5 一 一 c 
水 分 利用 RDI 初始 根深 10 cm m-c 
RRI 根深 每 日 最 大 增 量 1.2 cm/d c 
z= RDMCR 最 大 根深 120 cm m-c 
à 3.2 田间 验证 表 2 总 地 上 生物 量 和 时 面积 指数 验证 结果 
N 在 田间 试验 区 域 ， 通 过 2016 年 测量 D1 处 理 Table 2 Validated results of TAGP and LAI 
= 的 田间 数据 校正 模型 参数 ，2017 年 生长 季 内 ag TOPE E 
N D1~D5 处 理 的 时 间 序 列 的 总 地 上 生物 量 、 叶 面 。 处 理 e gna) NAMSE% RO RMSE NRMSEM 
o 积 指数 和 土壤 水 分 含量 被 用 于 验证 模拟 精度 。 T a ee a T 
B 3.2.1 TAGP 和 LAI 验 证 D2 097 874 11.1 0.95 0.17 10.9 
O TAGP 和 LAI 模拟 精 度 是 评价 束 树 生 长 模拟 D3 0.94 1084 16.9 0.94 0.15 10.0 
N ] - _ D4 0.92 1122 20.5 0.90 0.17 14.8 
in 过 程 的 重要 指标 。 表 2 显示 了 不 同 处 理 的 TAGP D5 0.92 869 19.8 0.79 0.18 21.1 
= 和 LAI 模 所 结果。 校正 模型 模 拟 的 D1~D3 处 理 enaa am MET 
X 的 TAGP 与 实测 值 一 致 性 较 好 (PCT 094), B we 
c 拟 精度 较 高 (8.7%< NRMSES16.9%), BEL D4 All Er E C DINE 
c D5 处 理 的 总 地 上 生物 量 的 性 能 低 于 Dl、D2 和 ”M5E 均 小 于 6.1%， 绝 对 误差 较 小 。 
o D3 A3, RH 0.92, NRMSE 约 等 于 20%。 模 拟 侍 合 总 地 上 生物 量 、 叶 面积 指数 和 土壤 水 分 


D1~D4 处理 LAI 的 一 致 性 和 精度 较 高 ，R 大 于 
0.9，NRMSE 小 于 15%，D5 处 理 的 模拟 结果 上 略 
35, RÆF 0.79, NRMSE 大 于 20%。 表 明 随 着 灌 
水 量 的 降低 ， 水 分 胁迫 强度 将 增加 ， 校 正 模型 对 
水 分 胁迫 响应 的 模拟 能 力 可 能 下 降 。 总 之 ,在 
70%~100% 对 照 灌 溉 量 时 ， 校 正 模型 都 表现 了 
较 好 的 总 地 上 生物 量 和 叶 面 积 指数 模拟 性 能 。 
3.2.2 土壤 水 分 含量 验证 

模拟 和 实测 D1~D5 处 理 的 土壤 水 分 含量 的 
散 点 图 见 图 2， 测 量 值 与 模拟 值 总 体 趋 势 保持 一 
致 ， 模 拟 与 实测 土壤 水 分 含量 的 尼 在 0.29~0.75 
之 间 ，RMSE 在 0.010—0.013 Z Ù], NRMSE 在 
4.1%~-6.1% 之 间 。D4 和 了 D5 处 理 的 模拟 和 实测 土 


含量 模拟 的 评价 指标 RR 和 NRMSE， 当 灌溉 量 不 
低 于 经 验 灌溉 量 的 80% 时 ， 模 型 具有 较 高 的 一 致 
性 和 模拟 精度 。 虽 然 在 70% 和 60% 的 经 验 灌 溉 
量 时 ， 模 型 模拟 束 树 生长 过 程 总 地 上 生物 量 和 叶 
面积 指数 的 性 能 略 低 ， 但 最 大 绝对 误差 
NRMSE<21.1%， 仍 可 以 实现 中 等 精度 的 模拟 。 
在 所 有 灌溉 处 理 中 ， 校 正 模型 都 表现 了 极 高 的 土 
壤 水 分 含量 模拟 精度 (NRMSE< 10% )。 


3.3 县 域 尺度 验证 


3.3.1 模拟 最 大 LAI 验 证 
分 别 输入 实测 55 个 吏 园 的 初始 束 树 干 重 参 
数 运 行 模型 ，55 个 束 园 的 最 大 LAI 模 拟 结果 散 点 
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图 2 D1~D5 处 理 的 土壤 水 分 含量 模拟 结果 


Fig. 2 Soil moisture simulation results for D1—DS treatments 


模拟 LAI 的 变异 系数 CF 值 分 别 等 于 21.6% 和 


图 见 图 3，2016 和 2017 年 模拟 最 大 LAI 与 实测 


LAI 的 尼 分 别 为 0.64 和 0.78， 
NRMSE 分 别 为 13.3% 和 10.7%, 


一 致 性 较 好 ， 
模拟 精度 较 高 。 
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(a)2016 年 模拟 和 测 


量 最 大 LAI 的 散 点 图 


21.2%， 表 达 了 显著 的 空间 变异 
始 束 树干 重 区 域 空间 分 布 变异 。 
2.6 7 
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(b)2017 年 模拟 和 测 


图 3 2016 和 2017 年 县 域 最 大 LAI 模 拟 结 果 


量 最 大 LAI 的 散 点 图 


Fig. 3 The simulated results of the maximum LAIat county scale in 2016 and 2017 


3.3.2 ”模拟 单产 验证 


2016 Fil 2017 年 模拟 与 实测 单产 的 一 致 性 指 


标 尼 分 别 为 0.48 和 0.60 (图 4)， 
RMSE 分 别 为 0.83 和 0.94 kg/ha, 


单产 预测 精度 
NRMSE 分 别 是 


12.1% 和 11.9%, 
(NRMSE<20% ) 。 


数 、 干 物质 分 配 、 


具有 和 较 高 


的 单产 预测 


性 ， 也 解释 了 初 


精度 


图 4 也 表明 高 产 区 域 明显 被 低 
估 ， 影 响 模 拟 产量 的 因素 较 多 ， 包 括 CO, 同化 参 


叶 面 积 与 干 重 比值 和 叶片 老化 
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速率 。 本 研究 结果 可 能 的 原因 是 叶片 老化 速率 参 
数 设置 值 过 低 ， 田 间 试 验 测 量 的 叶片 老化 速率 值 
在 40~60 之 间 ， 县 域 尺度 部 分 果园 的 叶片 老化 速 
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率 值 可 能 高 于 设置 的 固定 值 55， 甚至 可 能 高 于 
60， 所 以 导致 高 产 区 域 被 低估 。 
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(b)2017 模拟 和 测量 最 大 LAI 的 散 点 图 


图 4 2016 和 2017 年 县 域 单产 模拟 结果 


Fig.4 Yield simulation results at the county scale in 2016 and 2017 


单产 模拟 的 相对 误差 频率 分 布 如 图 5 所 示 ， 
大 部 分 样本 频率 分 布 在 y=0 附 近 ， 平 均 相 对 误差 
分 别 为 -3.23% 和 -2.94%， 相 对 误差 绝对 值 的 平 
均值 分 别 为 9.72% 和 9.68%， 单 产 预测 精度 较 高 。 
CV 分 别 是 12.5% 和 11.9%， 具 有 一 定 产量 空间 分 


相对 频率 /% 


-30 -20 -10 0 10 20 30 
模拟 单产 的 相对 误差 /% 
(a)2016 年 模拟 单产 的 相对 误差 分 布 
注 : Bin Size=2.5 ,曲线 平滑 为 B-spline 


布 差异 ,不 同 树龄 和 种 植 密度 导致 的 初始 京 树干 
重 的 空间 异 质 性 可 能 导致 产量 分 布 差异 。 然 而 ， 
2016 和 2017 年 分 别 有 61.8% 和 63.6% 的 样本 被 低 
估 ， 叶 片 老化 速率 空间 变异 性 可 能 是 导致 模拟 精 
度 下 降 的 原因 。 


相对 频率 /% 


-30 -20 -10 0 10 20 30 
模拟 单产 的 相对 误差 /% 
(b)2017 年 模拟 单产 的 相对 误差 分 布 


图 5 模拟 单产 的 相对 误差 频率 分 布 


Fig. 5 Relative error frequency distribution of the simulated yield 


3.3.8 ”模拟 实际 蒸 散 量 和 水 分 利用 效率 验证 
模拟 的 55 个 观测 果园 的 实际 葵 散 量 和 水 分 
利用 效率 与 测量 值 的 相对 百分比 误差 频率 分 布 如 
图 6 所 示 。 所 有 果园 的 实际 E7, 模 拟 误差 小 于 
20%， 其 中 49 个 果园 的 模拟 误差 小 于 10%。2016 


fll 2017 年 模拟 ET7, 均 方 根 误差 分 别 为 36.1 mm 
(7.9%) 和 30.8 mm (7.4%) ( 表 3)， 表 现 了 极 高 
的 模拟 精度 (NRMSE<10%). 2016 EREM ET, F 
均 相 对 误差 略 高 于 测量 值 ，2017 年 低 于 测量 值 ， 
平均 相对 误差 的 绝对 值 分 别 为 6.44% 和 5.9896, 
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相对 频率 /% 
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0 
模拟 ET, 的 相对 误差 /% 
(a) 模 拟 实际 蒸 散 的 相对 误差 分 布 


图 6 模拟 实际 蒸 散 和 水 分 利 
Fig. 6 Relative error frequency distribution of the simulated ET, and WUE 


校正 模型 的 整体 E7, 模 拟 精度 较 高 。 模 拟 的 WUE 
相对 误差 略 低 于 ET,， 在 2016 和 2017 年 两 个 生长 
季 ， 模 拟 WUE 的 相对 百分比 误差 小 于 10% 的 样 
本 分 别 占 50.9% 和 47.3%， 小 于 20% 的 分 别 占 
85.5% 和 80%， 最 大 值 分 别 为 37.4% 和 38.3%。 
2016 年 模拟 WUE 的 平均 相对 误差 低 于 测量 值 ， 
2017 年 略 高 于 测量 值 ， 平 均 相 对 误差 的 绝对 值 分 
别 为 11.39% 和 12.73%。 校 正 的 模型 表现 了 较 高 
AY WUE 模拟 精度 (10%<NRMSE<20%), RMSE 
值 分 别 为 0.23 和 0.28 kg/m’. 
RI 模拟 实际 蒸 散 和 水 分 利用 效率 的 验证 指标 


Table 3 Validated indices of the simulated actual evapotrans- 


piration and WUE 
绝对 误差 平均 
参数 Æ RMSE NRMSE/% MBE/% 
值 M4E/% 
2016 36.1 mm 7.9 0.68 6.44 
T, 
2017  30.8mm 7.4 -2.80 5.98 
2016 0.23 kg/m? 15.3 -3:23 11.39 
E 
2017 0.28 kg/m? 14.8 0.68 12.73 
ik: MBE 为 相对 误差 的 平均 值 
+, * 
4 tt it 


4.1 田间 尺度 模拟 精度 分 析 

田间 尺度 验证 结果 表明 ， 基 于 时 间 序 列 的 总 
地 上 生物 量 和 叶 面 积 指数 的 模拟 结果 均 达 到 中 等 
精度 以 上 ， 土 壤 水 分 含量 达到 极 高 的 模拟 精度 。 


i == -Line x-0Ü 
1 2016 年 

1 2017 年 
I 


相对 频率 /% 


模拟 WUE 的 相对 误差 /% 
(b) 模 拟 水 分 利用 效率 的 相对 误差 分 布 


用 效率 的 相对 误差 频率 分 布 


SRI, DI—D3 的 模拟 性 能 略 高 于 D4 和 了 D5 处 理 ， 
说 明 校正 的 模型 随 着 灌溉 量 的 降低 、 水 分 胁迫 的 
增强 ， 模 拟 精 度 可 能 呈现 下 降 的 趋势 。 原 因 可 能 
是 研究 中 仅仅 使 用 2016 年 D1 处 理 数据 校正 模型 
参数 ，D1 处 理 接 近 潜 在 水 分 供应 ，2017 年 
D1~D3 处 理 灌溉 量 大 于 80%CK, 与 2016 年 D1 
处 理 的 灌溉 量 差 异 不 大 ， 现 树 面 临 水 分 胁迫 较 
小 ， 对 现 树 产量 影响 也 较 小 。D4 和 D5 处 理 可 能 
导致 更 严重 的 干旱 胁迫 ， 实 际 生 产 中 对 束 树 生长 
的 影响 可 能 高 于 模拟 过 程 的 干旱 胁迫 修正 结 
进而 导致 总 地 上 生物 量 和 LAI 的 模拟 误差 。LAI 
的 模拟 偏差 也 将 导致 D4 和 D5 人 处理 的 束 树 蒸腾 和 
EAR AW, PET REL BOK ET 
拟 精 度 略 低 于 D1~D3 人 处 理 。 

另外 ， 在 大 田 试验 过 程 中 ， 束 树 生 长 通常 会 
收 到 施肥 、 病 虫害 、 风 、 热 胁迫 等 因素 影响 。 在 
本 研究 中 ， 模 型 无 法 响应 这 些 因素 对 束 树 生长 的 
影响 ， 在 已 校正 的 衷 树 生 长 模型 中 耦合 经 典 的 氮 
素 运 移 模 块 或 病虫害 响应 模型 有 望 进一步 提高 模 
拟 精 度 和 机 理性 。 


4.2 模型 区 域 尺 度 应 用 的 不 确定 性 


在 运行 WOFOST 作 物 生长 模型 进行 生长 模 
拟 和 产量 评 佑 时， 输入 参数 、 气 象 驱 动 数据 的 不 
确定 性 以 及 模型 模拟 过 程 的 简化 都 会 影响 单产 估 
算 的 准确 性 汪 。 本 研究 预期 模拟 误差 主要 来 自 
输入 参数 的 不 确定 性 ， 初 始 束 树干 重 参 数 强 烈 影 
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响 农作物 的 初始 生长 速率 ， 并 显示 出 高 度 的 不 确 
定性 ， 这 可 能 会 影响 初始 LAI 和 最 大 LAI 的 生长 
速率 5 ， 进 而 影响 生长 模拟 和 单产 评估 精度 。 
通过 实测 55 个 果园 的 初始 束 树 干 重 值 驱动 模型 
在 县 域 尺度 模拟 最 大 LAI 和 单产 精度 较 高 ， 说 明 
县 域 尺 度 的 束 园 初始 束 树 干 重 空间 异 质 性 较 大 。 
在 已 报道 的 研究 中 ， 当 使 用 相同 树龄 果园 的 初始 
枣 树干 重 平均 值 驱动 模型 时 ，181 个 区 域 束 园 的 
单产 验证 NRMSE 为 16.3% (2016 年 ) 和 17.2% 
(2017 年 )"”， 明 显 高 于 本 研究 中 的 12.1% 和 
11.9% 的 验证 误差 ,这 也 说 明 初 始 束 树干 重 为 束 
树 生 长 模拟 过 程 的 高 敏感 参数 。 男 外 ， 尽 管 决 定 
叶片 衰老 速率 和 时 间 的 叶片 老化 速率 参数 主要 取 
决 于 作物 遗传 学 特性 ， 但 该 参数 也 受 干旱 、 营 养 
胁迫 、 病 忠 害 以 及 作物 管理 等 因素 的 影响 。 
但 WOFOST 模 型 无 法 模拟 这 些 因素 对 叶片 老化 
速率 的 影响 站。 研究 中 为 叶片 老化 速率 设置 一 
个 固定 值 (50) 无 法 表达 所 有 55 个 样本 的 实际 
叶片 老化 或 损伤 情况 ， 模 拟 最 大 LAI 的 误差 最 可 
能 来 自 于 叶片 老化 速率 的 空间 分 布 误 差 。 在 以 往 
研究 中 ， 通 过 遥感 反 演 LAI 同 化 潜在 模式 的 WO- 
FOST 模 型 可 以 降低 初始 束 树 干 重 和 叶片 老化 速 
率 的 不 确定 性 ， 提 高 区 域 尺度 的 更 树 单产 模拟 精 
度 ， 集 合 卡尔 曼 滤 波 EnKF 和 SUBPLEX 算 法 同 
化 遥感 反 演 LAI 后 模拟 单产 的 NRMSE 分 别 为 
9.2% 和 8.3% ”。 表 明 初 始 束 树干 重 和 叶片 老化 
速率 为 区 域 尺度 吏 树 生长 模拟 的 不 确定 性 参数 ， 
遥感 同化 方法 可 以 被 考虑 降低 这 一 不 确定 性 。 
另外 ， 区 域 尺 度 的 土壤 属性 的 空间 异 质 性 也 
会 导致 最 大 LAI 和 水 分 运 移 过 程 的 模拟 精度 。 同 
化 遥感 反 演 的 LAI 和 ET 状态 变量 可 以 降低 土壤 
水 分 含量 模拟 的 不 确定 性 ， 进 而 提高 产量 评估 精 
度 。 因 此 ， 县 域 尺度 单产 和 莹 散 量 模拟 误差 可 能 
来 自 灌溉 、 施 肥 和 修剪 管理 导致 的 初始 现 树 干 
重 、 叶 片 老化 速率 和 土壤 属性 参数 的 不 确定 性 。 
在 今后 的 研究 中 ， 可 通过 遥感 同化 土壤 水 分 、 著 
发 蒸腾 和 叶 面 积 指数 等 状态 变量 对 区 域 尺度 的 初 
始 束 树干 重 、 叶 片 老化 速率 和 土壤 属性 参数 进行 


校正 ， 提 高 区 域 尺度 模拟 精度 ， 为 在 区 域 尺度 分 
析 水 分 胁迫 对 束 树 生长 和 产量 的 影响 提供 机 理性 
模拟 模型 。 


5 结论 


本 研究 以 WOFOST 模型 为 基础 ， 以 县 域 尺 
度 的 束 树 生长 模拟 和 水 分 利用 评价 为 目标 ,使 用 
田间 实测 束 树 生理 生化 指标 和 土壤 理化 数据 校正 
模型 参数 。 在 田间 试验 果园 获得 了 较 高 的 TAGP 
和 LAI 模 拟 精 度 (8.3% <NRMSE <21.1%)， 以 及 
极 高 的 土壤 水 分 含量 模拟 精度 
(4.1% <NRMSE<6.1%)。 使 用 校正 的 模型 和 55 个 
果园 测量 的 初始 吏 树 干 重 驱动 模型 进行 了 县 域 尺 
度 的 束 树 生长 模拟 和 水 分 使 用 评价 。 结 果 表 明 ， 
2016 和 2017 年 模拟 最 大 LAI 和 单产 的 NRMSE 均 
小 于 20%， 实 现 了 较 高 的 模拟 精度 。 校 正 模型 也 
表现 了 极 高 的 ET, 模 拟 精 度 (NRME<10%) 和 较 
高 的 WUE 模 拟 精度 (NRMSE<20%)。 总 之 ,在 
田间 和 区 域 尺 度 的 刺 园 ，WOFOST 模 型 都 表现 
了 较为 理想 的 更 树 生长 模拟 精度 和 水 分 利用 评价 
性 能 ， 可 为 分 析 和 气象、 土壤、 灌溉 管理 和 刺 树 生 
长 耦合 影响 提供 一 种 新 思路 。 
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Dynamic Simulation of Jujube Tree Growth and Water 
Use Evaluation Based on the Calibrated WOFOST Model 
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University, Beijing 100083, China) 


Abstract: Irrigation schemes determined based on statistical analysis of field trials are usually only applicable to specific soils 
and meteorological environments. It is difficult to quantitatively analyze the impact of irrigation strategies on the growth of ju- 
jube trees. In order to realize the quantitative analysis of the influence of temperature, light and water resources on the growth of 
fruit trees, WOrld FOod Studies (WOFOST) model parameters were calibrated to simulate the jujube tree growth and water mi- 
gration process. Firstly, the observed data obtained from field trials in 2016 and 2017 were used to calibrate the phenology de- 
velopment, initialization, green leaf, CO, assimilation, dry matter partitioning, respiration, and water use parameters of the WO- 
FOST model. Secondly, the time series of total above-ground biomass, leaf area index (LAI) and soil moisture content in field 
trials were dynamically simulated, and accuracy verification and analysis were also performed. Finally, the maximum LAI, 
yield, actual evapotranspiration (ET, ) and water use efficiency (WUE) data of 55 orchards were employed to evaluate the per- 
formance of the calibrated model at the county scale. The results showed that the coefficient of determination R’ of TAGP simu- 
lated in the field test area was between 0.92 and 0.98, and the normalized root mean square error (NRMSE) was between 8.796 
and 20.5%, the R of simulated LAI ranged from 0.79 to 0.97, and the NRMSE ranged from 8.3% to 21.1%. The R of the simu- 
lated soil moisture content was between 0.29 and 0.75, and the NRMSE ranged from 4.1% and 6.1%. The model could well sim- 
ulate the time series of jujube tree growth dynamics and soil moisture content changes. At the county scale, the R^ between the 
simulated and measured maximum LAI were 0.64 and 0.78, and the NRMSE were 13.3% and 10.7% in 2016 and 2017, respec- 
tively. The simulated yield showed R° value of 0.48 and 0.60, and NRMSE of 12.1% and 11.9%, respectively. RMSE of the 
simulated versus measured ET, were 36.1 mm (7.9%) and 30.8 mm (7.4%), respectively. The model also showed high WUE 
simulation accuracy (10%<NRMSE<20%) with RMSE values of 0.23 and 0.28 kg/m in 2016 and 2017, respectively. In short, 
WOFOST model achieved accurate simulation of jujube tree growth and water transport at the field and county scales, which 
may provide new ideas for the quantitative and mechanism analysis of the coupled effects of soil, weather, irrigation manage- 
ment and jujube tree growth. 


Key words: crop growth model; parameter calibration; WOFOST; jujube tree; water use efficiency 
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